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We discuss the spin relaxation of a strongly correlated two-dimensional 
(2D) electron gas (2DEG) in the quantized Hall regime when the filling factor 
is close to an odd-integer. As the initial state we consider a coherent deviation 
of the spin system from the B direction and investigate a break-down of this 
Goldstone-mode state due to the spin-orbit (SO) coupling and smooth disor- 
der. The spin relaxation (SR) process is considered in terms of annihilation 
transitions in the system of spin excitons (magnons). 
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Great bulk of recent SR measurements and theoretical studies deal with 2D electrons 
confined in quantum wells or dots. 1 10 However, during the last two decades only a few of 
works have been devoted to the SR in the quantized Hall regime proper. 1-6 Specifically, we 
imply conditions under which a 2DEG is in a strong perpendicular magnetic field (B > 10 T) 
and in the absence of holes and magnetic impurities. We exclude from the consideration the 
effects of electrons in the edge states (cf. Ref. 6), and restrict our study to an odd-integer 
filling. 

In the studied problem the relevant SR time is actually not a spin dephasing time but 
a time of Zeeman energy relaxation conditioned by a spin-flip process. Indeed, any spin-flip 
means actually a reduction of the Zeeman energy \gfi B BAS z \ (B \\ z, AS Z = S z — S is the 
S z component deviation from the equilibrium value So, g ~ —0.44). The mechanism, which 
makes the relaxation irreversible, has thereby to provide the energy dissipation. Another 
necessary condition is a spin-flip mechanism non-conserving S z . 

In the present work the SO coupling is considered as the cause mixing spin states and the 
disorder [to be more precise, the smooth random potential (SRP)] as a dissipation. A crystal 
lattice is still implicitly assumed to be present as a "cooler" for 2DEG excitations. We will 
suppose that all electron-phonon relaxation processes which do not change the 2DEG spin 
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state occur much faster than the SR. [Thermodynamic relaxation times are estimated to be 
<lns, i.e. well shorter than the SR time (see the corresponding estimate given at the end 
of the paper).] 

We thus study the case where the 2DEG is a quantum Hall "ferromagnet" (QHF); i.e. 
the filling factor is v — N/N^ ~ 2k + 1, where Af and = L 2 /2nl 2 B are the numbers of 
electrons and magnetic flux quanta [L 2 is the 2DEG area, Ib is the magnetic length). In the 
high magnetic field limit, which really represents the solution to the first order in the ratio 
r c = (e 2 /elB)/hu c considered to be small (u c is the cyclotron frequency, e is the dielectric 
constant), we get the ground state with zeroth, first, second,... and (k— l)-th Landau levels 
fully occupied and with nth level filled only by spin-up electrons aligned along B. Only a 
bare handful of experimental results on the SR in such a QHF were obtained, first indirectly 
(the line-widths of the electron spin resonance (ESR) were measured in Ref. 1) and then 
directly (in communication of Ref. 2 the photoluminescence dynamics of spin-up and spin- 
down states was studied). The measured times are 5 — 10 ns and exceed by 1-2 orders the SR 
times observed in quantum wells where the SR process is governed by the spin interaction 
with band holes. 7 

As the spin-system perturbation in the QHF we study a coherent deviation, when the 
total S number is not changed. Namely, the initial state is a Goldstone mode which represents 
a quantum precession of the vector S around the B direction: \i) = ( S- ] |0). Here |0) stands 
for the QHF ground state and 5*- = V. is the lowering spin operator [j labels electrons; 
<j± = (& x ±ia y )/2, where cr x ,y,z are the Pauli matrices]. The number N is assumed to be 
macroscopically large: O^NkN^. The spin numbers of the \i) state are S=S =N ( p/2 and 
S z =N^/2-N (i.e. AS z \ t=Q = -N). The total Hamiltonian has the form H tot ^J2j H i ] '+H hlt . 
Here ifi nt is the many-electron (Coulomb interaction) part of the Hamiltonian which has the 
usual form (see, e.g., Ref. 11), and H 1 is the single-electron part: 



where q = — iV + eA/ch is a 2D operator, ez = \g\fisB is the Zeeman energy of one spin- 
flipped electron, and ip(r) is the SRP field [r has components (x,y)}. The SO Hamiltonian 
is specified for the (001) GaAs plane, 



#i = ^ 2 q 2 /2m: - e z a z /2 + H so + tp(r) , 



(1) 



H so = ol (q x cr) z +!3 {q y cr y -q x a x ) , 



(2) 
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and presents a combination of the Rashba term 12 (with the coefficient a) and the crystalline 
anisotropy term 13 (see also Refs. 3-5, 8). The parameters a and (3 are small: a, (3 <C IbTujJc 
(moreover, really a</3~10 -7 K • cm<l B ez)- This enable us to account H so perturbatively. 

If the SRP is assumed to be Gaussian, then it is determined by the correlator K(r) = 
(tp(r)tp(0)) , where tp(r) is the SRP field. We choose also (y?(r)) = 0, i.e. the SRP energy 
is measured from the center of the Landau level. In terms of the correlation length A 
and Landau level width A, the correlator is K(r) = A 2 exp (— r 2 / Is 2 ). In the realistic case 
A Pa 5 - 10 K, and A ~ 30 - 50 nm. We will study the case 

T < T* < e z < A < e 2 /el B < hu c , and A > l B (3) 

(T is the temperature which is actually assumed to be zero in the calculations; the value T* 
will be defined subsequently). All results which follow are obtained in the leading approxi- 
mation corresponding to these inequalities. 

The QHF is also remarkable for the following reason: when neglecting the last two 
terms in Eq. (1) the spin excitons (SEs) are actually exact (to the first order in r c ) lowest- 
energy eigen states. The most adequate description of the SE states is realized by the SE 
creation operators. 3 ~ 5 > 14 However, before writing them out we choose the bare single-electron 
representation. As previously 4 ' 5 we use the spinor basis which diagonalizes the first three 
terms in the Hamiltonian (1) to the first order in u = (3\/2/IbJiuj c and v = a\f2/lBhuo c : 



\vVK+llfJ K+lp + lU^/Klp K - lp J V Wk P J 

where i\) K , p = L~ l / 2 e wy (t) K {pl 2 B + x) is the wave function of an electron in the Landau gauge 
((f> K is the harmonic oscillatory function). We note that here and in the following we present 
only the perturbation expansion to within the framework of the leading order in u and v. 
The exciton creation operator is 



where a p and b p are the Fermi annihilation operators corresponding to the states (4). The 
annihilation excitonic operator is Q abq = ^<2jj 6q j = Q^-q an< ^ we em pl°y a l so the "shift" 
operators A{ = iV^ 1/2 QL q and B{ = Nf^Q^ (A n = A^, £ q = B^.) In Eq. (5) the 
orbital index k is dropped since for our purposes the approximation of projection onto a 
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single Landau level is quite sufficient. In the following we drop also the "spin-orbit" index ab 
at the operator (5). Such single-level exciton operators constitute a Lie sub-algebra which is 
a part of an Excitonic Representation (ER) algebra (e.g., see Ref. 15 and references therein). 
In our case the relevant commutation rules are as follows: 

[Qqi> 2q 2 ] = A™ -e-^B qi ^ 2 , and 

(6) 

where 6 12 = ^(qi x Q2)z/2- Besides evidently [<2 qi , <2 q2 ] = [^4 qi , £> q2 ] = 0. The ground state 
|0) is completely determined by the equations AJO) =5 q ,o|0) and £> q |0) = 0. Single-exciton 
states are normalized: (0| <2 qi <2 q2 10) = <5 qi , q2 . 

In the limit A — > 0, - ► and at r c 1 the state 

|iV;l;q) = Q q (Qt^|0) (7) 

is the eigen state of the system studied. It has the spin numbers S = N^/2 — l and S z = 
N^/2—l—N [see below the expressions (8) which should be used to calculate S and S z ]. The 
corresponding energy is (N+l)ez + £ q , where £ q is the exchange part of the SE energy. The 
small momentum approximation ql B <1 is quite sufficient for our problem, and therefore 
£ q =(ql B ) 2 /2M^ K (general expressions for the 2D magneto-excitons can be found in Ref. 11). 
Here M XjK is the SE mass at v = 2k+1, namely: l/M Xj0 = (e 2 / el B )y r n~/%, 1/M x> i = 7/4M X>0> ... 
Spin operators in terms of the ER are invariant with respect to Hso'- 

S z = N <p (A -B )/2, S- = N^ 2 Ql S 2 = A^QjSo + S 2 + S z , (8) 

At the same time in the basis (4) the operators tp (r) and H int acquire corrections proportional 
to u and v. Specifically, calculating J ^(p(r)^d 2 r, where \I> = J2 p (a p ^ Kpa + b p ^ Kpb ), we get 
the terms responsible for a spin-flip process: 

<p = N l J 2 l B ¥>(q) (iuq+ - vq-) Q q + H.c. (at ql B < 1) . (9) 
q 

Here </?(q) is the Fourier component [i.e. <p = X] q ^(q) e * qr ]) an d <?± = = R(<2 , x ±iq y )/V2- 

We stress one essential feature of the states (7). In spite of the existence of a formal 
operator equivalence Qq = lim we find that \N;1;0) and lim \ N; 1; q) present different 

q— >0 H q^O 

states. Indeed, in these states the system has different total spin numbers S = N^/2 and 
S = N^/2 — l respectively. [S z = N^/2 — 1 — iV is the same for both.] So, the excitation of a 
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"zero" exciton (with zero 2D momentum) corresponds to the transition S z — > S z — 1 without 
any change of the total S number, but each "nonzero" SE changes the spin numbers by 1: 
S—*S — 1, S Z ^S Z — 1. Let us introduce the notation \N) = (So) |0). The initial state \i) 
in the ER is actually a "Goldstone condensate" (GC) containing N zero excitons: |i) = |iV). 

Our goal is to study the process of the GC break-down. The state \N) is certainly 
degenerate and we will solve the problem in terms of the quantum-system transitions within 
a continuous spectrum. The transition probability is determined by the Fermi Golden Rule: 
Wfi = (2n/h)\M.fi\ 2 5(Ef — Ei). In our case evidently the final state |/) is obviously the state 
where a part of the Zeeman energy has been converted into an exchange energy. Since it is 
exactly the single-electron terms that constitute the perturbation responsible for the M. fi 
matrix element, we could find that such a transition is the 2Xq— >X q * process in the lowmost 
order of the perturbative approach (we denote the zero exciton by X and the nonzero one 
by Xq). In other words the final state for this transition is |/) = \N — 2; 1; q*). The value q* 
is determined by the energy conservation equation Ef — Ei which reads 2ez = ez+£(q*), i.e. 
q* = y/2M^ K e z /l B . 

We point out that the SO interaction (2) alone does not provide a quantum fluctuation 
from the GC to any state with a different electron density (in particular, to the state \N — 
2; 1; q*)). This feature may be verified by means of direct analysis of the SO influence on the 
2DEG spectrum (c.f. Ref. 4). Besides, this can be recognized from general considerations. 
Indeed, each of the components % (i = x, y) commutes with any component of the operator 
P — ^2j[hcij — (e/c)B x r.,]. The latter commutes with the Hamiltonian ft 2 q 2 /2m* + i/i n t 
and one can also check that P, Q^j =qQ\ [if Q\ is defined by Eq. (5) to within the zero 
order in H S o\- The operator P hence plays the role of 2D momentum in the magnetic field. 16 
Its permutability with the total 2DEG Hamiltonian reflects the spatial homogeneity of the 
system under consideration in the "clean limit". Since [H so ,P] = 0, the SO coupling by 
itself cannot destroy the homogeneity in any order of the perturbation approach. 

The transition is thereby determined by the operator (9). The corresponding ma- 
trix element M if = ((q*; 1; N -2\<p\\N))[R(N)R(N -2; 1; q*)]~ 1/2 was actually calculated 
in Ref. 4. (The SRP plays the same role as the phonon field studied there.) Here and 
in the following the notation R(...) stands for the norm of the state !•••)• The result 
is | .Mi/ 1 2 = JV ^" 1 - > (u 2 + v 2 ) |<?*/M?(q*)| 2 , and we obtain the rate of the % — > / transition 
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(2n/k)Z (l \M tf \ 2 5(q 2 l 2 B /2M x , K -e z )=N(N-l)/TN^ (N > 1), where 

1/r = 8n\a 2 +^)Ml K e z K(q*)/h 3 u; 2 c l% . (10) 

Here K stands for the Fourier component of the correlator: K(q) = L 2 |^(g)| 2 /47r 2 . 

The quantum transition to the state |/) = |7V — 2; 1, q*) is certainly a first step in the SR 
process. This state is thermodynamically unstable. In a time which is much shorter than r 
it turns to a state |iV — 2; l,qo), where q takes the lowest possible nonzero value. In fact, 
relevant values of go are determined by the SRP field. The SE interaction with the SRP 
incorporates the energy C/ X . SRP ~g/^A/A (the "nonzero" SE possesses the dipole momentum 
elg[qxz], see Ref. 11). The latter determines the inhomogeneous uncertainty of the SE 
momentum 5q ~ M X)K A/A (this follows from the equation dqdS^/dq = f/ x -sR.p)- Therefore 
"quasi- zero" wave numbers in the range defined by inequalities 

0<g <M XiK A/A, (11) 

present the lowest physical limit for momenta of the nonzero SEs. Under the conditions (3) 
we find that q <C q*, and the SE exchange energy (q l B ) 2 /2M XtK is smaller than the value 
T*=M X)K (AIb/A.) 2 . This in turn is negligible in comparison with ez- 

To solve the problem in a complete form we have obviously to study the general state of 
the type: 

\N; Ml ,M 2 ,...,M K ) = (QL) Ml (2lJ M2 - (QU^)- (12) 
All the wave- vectors q ok are assumed to satisfy the condition (11). We will also use for this 
state a shorthand notation \N;M), where M = Y2k ^ k * s ^ ne total number of the nonzero 
SEs. If M 3>1, we assume that lCi^C N^. In the framework of our approach the state 
(12) is an approximate eigen state of a QHF having energy (N + M)ez and spin numbers 
S z = N ( f ) /2—N—M and S = N ( p/2—M. This value of S z is the exact one. It can be calculated 
employing the representation of Eqs. (8) and the commutation rules (6). The same algebra 
allows us to find that S|N;M) = S(S+1) (\N; M) + \s)), where the norm of the state \e) 
is small compared with the norm of \N;M), namely: R(e)/R(N; M) = 0(m 3 n/ K). The 
notations n — N/N^ and m — M/N^ are used for "reduced" quantum numbers. In the special 
case when N = 0, the state |0, M) can be treated as a "thermodynamic condensate" (TDC) 
which arises if M is larger than the critical number of nonzero SEs. The latter is estimated 
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at iV^/ f^/jexp [(£q + \U x -sKp\) /T] — 1} (e.g., c.f. Ref. 4), and in our case (3) it is at least 
smaller than N^M^T. At the same time, M is determined by the spin S of the system, 
therefore at a given M=N <t> /2—S>N ( j ) M XtK T we find that below some threshold temperature 
the nonzero SEs necessarily form a TDC. For macroscopically large iV and M the state (12) 
hence features a coexistence of GC and TDC. It should also be noted that specific values qofc 
as well as specific distribution given by M k numbers have no physical meaning. The final 
results should not depend on them but only on M and N. 

We can now write the kinetic equations corresponding to the relevant spin transitions. 
The rate dN/dt is determined by the 2X — > X q . (— > X qo ) process (which presents a GC 
depletion with a simultaneous "flow" to TDC) and by the X + X qo — > -X" q *(— > Xq^) one. 
The rate dM/dt is also formed by the 2X — >JT q *(— >X qo ) transition (which provides a TDC 
evolution) and by the X qo +X q ^— >JT q *(— >JT q ») one (determining a TDC depletion). [Values 
of qo,q' Q and belong to the region (11).] The corresponding equations are derived again 
with help of the Fermi Golden Rule and Eqs. (6) (with vanishing 9 12 in the latter): 

dn/dt = -(2fi nn + /w)/r, and dm/dt = (p nn - ^ mm )/T , (13) 



where 

rnn 



= E 



|(M;iV-2|Q q »Q^|iV; M)| 2 _ N A R{N-2; M+l; q*) 

i2(7V;M)i2(JV-2;M+l;q*) ~ N^R(N; M) 
\(M 1 ,...M k -l,...M K ;N-l\Q^ ok Q_ <l ,\N;M)\ 2 



1+0 



/,'S.Y: M)R(N — 1; M; q*) 
4^ W -l;M;q-_) [1+0 



7V|i?(7V; M) 

V- I (Mi, ...M fc -1, ...Mi-1, ...M K ;7V|Q Q^|iV;M)| 2 
/W ^ i?(iV; M)R(N] M-l; q*) 

^ 2.U'/,'(.V:.U-1: (| ) 

= 7V 2 i?(iV; M) [1+ ° (Jf/JV * )] 

[i2(iV; M + l; q*) is the norm of the <2 q *|iV; M) state]. In this way we find that the norms 
appearing in these equations satisfy the conditions R(N;M + l;q*)/R(N;M) = r, R(N + 
1; M)/R{N; M) = Nr n and R(N; M±, ...M fc +1, ...M K )/R{N; M) = M fe r m , so that 

/inn = n 2 r/rl, [i nm = Amnr/r n r m , fx mm = 2m 2 r/r 2 n , (14) 

where the factors r, r n , r m are determined by the equations 

1 = (1 - n - 2m) /r n + 0(m 3 /nK), 1 = (1 - 2n - 2m)/r rn + n 2 /r 2 + 0(m 2 ), 

(15) 

1 = (1 - 2n - 2m) /r + Amn/r n r m + n 2 /r 2 + 2m 2 /r 2 m + O(l/iV ) . 



(Positivity of r n pa l—n — 2m provides the physically obvious requirement \S Z \ < S.) The 
last terms in Eqs. (15) would just depend on the specific set of the numbers. We 
can therefore calculate r m and obtain the final result only in the m <C 1 case. Meanwhile, 
the values of m(t) are determined by the initial value n(0). According to Eqs. (13)-(15) 
max(m) ~n(f )[1 — n(t*)]/y/2 (t* is the time at which m peaks), i.e. at least m 2 <l/32. [In 
particular, at n(0)=0.5 max(m)~0.1.] The problem (13)-(15) is thus solved to the leading 
order in m. We should put r=r n =r m =l in /i mn and fi mm , but r n =l— n and r=(l— n) 2 in /i nn . 
This yields the analytical result n(t) = l/[2n(0)(t/r) 2 + 2(t/r) + l/n(0)] and m = n(t)n(0)t/r. 

The dependences are shown in Fig. 1. The vector S(t) at moments t = 0, r, 2r, ... 
is depicted in the inset. The time (10) hence governs the breakdown of the GC. (Under 
the realistic conditions above, r ~ 10 -8 — 10~ 7 s.) However the SR occurs certainly non- 
exponentially and the actual time is increased by a factor of N^/ AS z (0). 

Experimentally, the GC could probably be created by microwave pumping at the electron 
frequency ez/fr- This should cause to "rotate" the QHF spin without changing of the S 
modulus. As to observing of the SR variation with time, one can think that the optical 
technique 2 is relevant in the case. 
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FIG. 1. Time dependences of \&.S Z \/N ( f ) =n{t}\m{t) and of \/S.S\/N ( f,=m{t) are shown in the main 
picture for n(0)=| A52(0)|/^=0.455. The vectors S(i) at equidistant moments of time are plotted 
in the inset with step r. The dotted line is the arc with radius So=N ( p/2. The gap between the 
dashed and dotted lines reflects the deviation of the spin modulus from the value of Sq. 
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